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Abstract: Recent SNO results strongly favour the large mixing angle (LMA) MSW solar 
solution. We argue that there are only two technically natural low energy neutrino mass 
matrix structures consistent with the LMA MSW solution, corresponding to either a hier- 
archy or an inverted hierarchy with pseudo-Dirac neutrinos. We construct the MNS matrix 
to leading order in the small angle ^13 including the neutrino and charged lepton mixing 
angles and phases, the latter playing a crucial role for allowing the inverted hierarchy case 
to be consistent with the LMA MSW solution. We then consider the see-saw mechanism 
with right-handed neutrino dominance and show how the successful neutrino mass matrix 
structures may be constructed with no tuning and with small radiative corrections, leading 
to a full, partial or inverted neutrino mass hierarchy. In each case we derive approximate 
analytic relations between the input see-saw parameters and the resulting neutrino masses, 
mixing angles and phases, which will provide a useful guide for unified model building. 
For the hierarchical cases the LMA MSW solution gives a soft lower bound {Uesl ~ O-^i 
just below the current CHOOZ limit. Both hierarchical and inverted hierarchical cases 
predict small PPou with \mee\ ~ 0.007 eV within the sensitivity of future proposals such as 
GENIUS. Successful leptogenesis is possible if the dominant right-handed neutrino is the 
heaviest one, but the leptogenesis phase is unrelated to the MNS phases. 
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1. Introduction 

Recent SNO results on the neutral current (NC) flux and the day-night effects Q, when 
combined with other solar neutrino data especially that of Super-Kamiokande Q strongly 
favour the large mixing angle (LMA) MSW solar solution [Q] with three active light neutrino 
states, and 9i2 ~ 7r/6, Arn^i ~ 5 x 10~^eV^ |5|, [^. The atmospheric neutrino data is 
consistent with maximal — u-r neutrino mixing with |Am32| ~ 2.5 x 10~^eV^ and the 
sign of Am|2 undetermined. The CHOOZ experiment limits ^13 ^ 0.2 over the favoured 
atmospheric range Q. The combined neutrino data is well described by an MNS matrix 
with 6^23 7r/4, 6^12 ^ ir/6, 6*13 < 0.2, which we refer to as the LMA MNS matrix. 
It is clear that neutrino oscillations, which only depend on Am?^ = mf — gives no 
information about the absolute value of the neutrino mass eigenvalues m-j. Recent results 
from the 2df galaxy redshift survey indicate that ^ < 1.8eV(95%C.L.) under certain 
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mild assumptions |11|. Combined with the solar and atmospheric oscillation data this 
brackets the heaviest neutrino mass to be in the approximate range 0.04-0.6 eV. The fact 
that the mass of the heaviest neutrino is known to within an order of magnitude represents 
remarkable progress in neutrino physics over recent years. 

The basic possible patterns of neutrino mass consistent with this data are: (i) hierarchy 
(full mi <^ 7712 1TT-3, or partial mi ^ m2 ^ w-s), (ii) inverted hierarchy (mi ~ m2 ^ m^, 
or the pseudo-Dirac form —mi ~ m2 ^ ^3)1 or (iii) degenerate m^ « m^ ~ m^ |12|. 



Although oscillation data does not distinguish between these possibilities, the theoretical 
requirement that the neutrino spectrum is generated in a technically natural way, does 
provide an additional guiding principle. For example it is clear that a degenerate mass 
scale m? ^ |Am?j| implies very small fractional neutrino mass splittings. For example, 
if ms = O.SOOeV then atmospheric oscillations require neutrino masses m2 = 0.497eV. 
The problem is one of tuning, both to set up the small mass splitting at the high energy 



scale, and to preserve it in the presence of radiative corrections |13|. In general such a 
neutrino spectrum is not technically natural, since small perturbations in the high energy 
input parameters will violate the low energy degeneracy. One can envisage a technically 
natural mechanism which would lead to a degenerate pair of neutrinos with opposite sign 
masses mi ~ — m2, but to achieve the full three- fold degeneracy is much more difficult 
|14],|T^. A similar objection can be raised against the inverted hierarchical spectrum in 
which the almost degenerate neutrinos have the same sign masses mi ~ m2 where for 
mi = O.OSOOeV we require m2 = 0.0497eV for the LMA MSW solution. In these kinds of 
inverted hierarchy models there is no natural mechanism that can set and preserve such 
mass splittings, however as we shall see, in the inverted hierarchy model with a pseudo- 
Dirac neutrino pair corresponding to opposite sign masses mi « — m2, a natural mechanism 
is possible. 

Assuming the LMA MNS matrix plus naturalness implies that there are then just 
two leading order forms for the light physical effective Major ana neutrino mass matrix 
^'ll (where LL means that it couples left-handed neutrinos to left-handed neutrinos) 
corresponding to either a hierarchical spectrum with m ~ ms ^ mi,m2 arising from: 



<L « 1 1 - (1.1) 




or an inverted hierarchical spectrum with a pseudo-Dirac pair of neutrinos with m ~ m2 
—mi ^ m3 arising from: 

/ 1 -T 



m 



LL 



10 0-= (1.2) 
-10 ' ^2 



where in both cases m ~ Y^|Am|^ ~ 5 x lO^^eV. It is remarkable that only by assuming 
the LMA MNS matrix plus naturalness arguments we are led to only two possible leading 
order structures for the neutrino mass matrix, both of which have a zero in the 11 element 
corresponding to a small /3/3oi/ rate corresponding to \mee\ = at leading order. 
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For the hierarchical case the experimental requirement |Am|2l ^ l^^^^iil implies that 
m| 3> m| which again looks technically unnatural since we would expect two roughly equal 
eigenvalues once the lower block is diagonalised. For the LMA MSW solution we only 
require a mild hierarchy, |?Ti2|/|m3| ~ 0.1, and this will require accidental cancellations 
of order 10% to take place in the diagonalisation of the lower block of the neutrino mass 
matrix. Although the degree of tuning necessary to achieve the hierarchy is not large, it 
does imply that the radiative corrections of the matrix elements will be competetive with 
the amount by which they need to be tuned, and hence that radiative corrections will be 
very important in determining the low energy spectrum. In the case of see-saw models 
1 16] the radiative corrections may be sufficient to destroy (or create) the cancellations 



necessary to achieve the desired hierarchy |17|, [p^ . ^ In all cases radiative corrections 
will be important and the neutrino masses and mixings calculated in a given high energy 
unified theory will not be simply related to the low energy ones. 

The above situation could be improved if it were possible to make the neutrino mass 
hierarchy completely natural. If the hierarchy ni^ ^ m2 could emerge without any tuning 
at all, not even at the level of 10% accidental cancellations, then the low energy spectrum 
would faithfully preserve the nature of the spectrum calculated at high energy, without 
being severely affected by radiative corrections. If this could be achieved then the low 
energy measurements would provide a direct window into the nature of the high energy 
theory. Therefore it is interesting to ask whether it is possible for the neutrino mass 
hierarchy to arise in a completely natural way? Indeed this is possible if the see-saw 
mechanism is supplemented by a mechanism known as single right-handed neutrino 



dominance [^], |23|, |24]. According to single right-handed neutrino dominance one of the 
right-handed neutrinos makes the dominant contribution to the lower block of m'^j^ causing 
its determinant to approximately vanish, and thereby leading to |m2| <C jmal without 
relying on accidental cancellations which are subject to important radiative corrections. 
Single right-handed neutrino dominance does not mean that there is only a single right- 
handed neutrino, only that one of the right-handed neutrinos is making the dominant 
contribution. If the dominant right-handed neutrino is denoted i'rs with a heavy Majorana 
mass Y and Dirac couplings d,e,f to the weak eigenstate neutrinos v^.-.i'ij.-.i't given by 
^Rz{d^e + ei^/i + fT^r) then according to the see-saw mechanism this will result in a light 
physical neutrino ~ dVf, + eu^ + fur of mass ms ~ (d^ + + f'^)/Y [22| together with 
two light orthogonal combinations of neutrinos which would be massless in the limit that 
there is only a single right-handed neutrino. The requirements of a maximal atmospheric 
angle 823 ~ vr/4 and a small CHOOZ angle 613 <^ 1 imply the relation d <^ e ^ f p^ . 

In order to account for the solar data we must consider the effect of the sub-leading 
right-handed neutrinos. These perturb the spectrum leading to a small second neutrino 
mass m| <C m^. The strength of the hierarchy is controlled by the relative importance of 
the sub-leading right-handed neutrinos, rather than relying on accidental cancellations. If 
the sub-leading contribution is dominated by a single sub-leading right-handed neutrino 
i^R2 with mass X and couplings i'R2{ai'e + bv^ + cvr)-, then this leads to a second neutrino 

^The one- loop beta function coefficients in the Standard Model have recently been corrected in |^^, [ plf . 
The two loop result for the MSSM is given in p^. 
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mass of order m2 ~ (l/2)(6 — cf' jX which only depends on the subleading parameters 
because we must have m2 = in the hmit of there being only a single right-handed 
neutrino [24|. The solar angle was given by tan0i2 ~ ^/2a/{h — c) Note that the solar 
angle is completely determined by the sub-leading couplings, due to a natural cancellation 
of the leading contributions [24|. The lightest neutrino mass mi is generated by the sub- 
sub-leading couplings due to the right-handed neutrino vri with mass X' leading to a 
full neutrino mass hierarchy mi ^ m2 <C m-3 [24|. We shall refer to this as sequential 
sub-dominance . 

The sub-leading contribution may alternatively result from two equally contributing 
sub-dominant right-handed neutrinos in which case only a partial neutrino mass hierarchy 
results mi 1z m2 <^ m^ and the results above will be different [24|. While the full hierarchy 
results from an approximately diagonal right-handed neutrino mass matrix, a partial hier- 
archy may result from three possible textures for the right-handed neutrino mass matrix 
namely diagonal, democratic or off-diagonal where the nomenclature refers to the upper 
block IH: 

/ X' 
0X0 

^ y 

( X X q\ 
X X 

Vo y 

/ X 

Mf-^'''^ = X 

Y y 

where the heavy Majorana neutrino mass matrices Mnn couple right-handed neutrinos 
to right-handed neutrinos. Note that though the right-handed neutrinos in Eq.L5 have 
a common pseudo-Dirac mass X, this may or may not lead to the lighter two physical 
neutrinos arising from the matrix in Eq.Ll having a pseudo-Dirac mass. As discussed [24| 
Eq jl.3| can lead to either a full or a partial hierarchy, while Eqs.1.4 and 1.5 only lead to 
partial hierarchies. We shall refer to the cases where partial hierarchies result from the 
structures in Eqs.L^, 1.4 and 1.5 as diagonal sub-dominance, democratic sub-dominance 
and off-diagonal sub-dominance. We shall continue to refer the case where a full hierarchy 
results from the diagonal structure in Eq.1.3 as sequential sub-dominance. Clearly they are 
all sub-classes of single right-handed neutrino dominance since in all cases the dominant 
contribution to the neutrino mass matrix comes from the single right-handed neutrino with 
mass y. The radiative corrections in such an approach would be expected to be small since 
the hierarchy is now completely natural, and this has been verified explicitly by a study 
of each of the cases in Eqs. [1.3|jL^J1.5| the result of which shows that radiative corrections 
typically change the physical predictions by only a few per cent [p5|. 



RR 



(1.3) 



(1.4) 



(1.5) 



The neutrino mass matrix in Eq.L2 corresponding to an inverted hierarchy with oppo- 
site sign masses mi ~ —772.2 and —mi ~ m2 ^ 772-3, can be reproduced by three right-handed 
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neutrinos with the texture in Eq.L5 [26|. However an important difference is that now it 
is the off-diagonal pair of right-handed neutrinos with pseudo-Dirac mass X that domi- 
nates the neutrino mass matrix leading to the pseudo-Dirac structure of neutrino masses 
in Eq jl.2| , with the right-handed neutrino of mass Y now giving the sub-dominant contri- 
butions [^6|. We shall refer to this as off-diagonal right-handed neutrino dominance. As 
in the hierarchical cases based on single right-handed neutrino dominance, the resulting 
inverted hierarchical spectrum does not rely on any accidental cancellations and is tech- 
nically natural. The radiative corrections arising to the inverted neutrino mass spectrum 
arising from off-diagonal right-handed neutrino dominance have been studied and shown 
to be only a few per cent [26|. 

In this paper we consider neutrino mass matrices with the leading order structures in 
Eqs.l.l,L^, and diagonalise each of them to leading order in ^13 to extract the neutrino 
masses, mixing angles and phases. The MNS matrix is then constructed to leading order 
in the small angle ^13 including the neutrino and charged lepton mixing angles and phases, 
the latter playing a crucial role for allowing the inverted hierarchy solution to be consistent 
with the LMA MSW solution. We then go on to show how the neutrino mass matrix 
structures in Eqs. |l.l|jL^ may be constructed naturally from the see-saw mechanism with 
right-handed neutrino dominance, with no tuning and with small radiative corrections, 
leading to a full, partial or inverted neutrino mass hierarchy. In each case we derive 
approximate analytic relations between the input see-saw mass matrices and the resulting 
neutrino masses, mixing angles and phases. The goal of this analysis is to provide a useful 
and reliable guide for constructing the LMA MNS matrix in unified models ||27||. 



The analysis builds on that in |24] and by including the effects of phases and 
the charged lepton mixing angles, which we did not previously consider. It is sufficient 
in a top-down approach to work to order ^13, since although the radiative corrections 
are only a few per cent (due to right-handed neutrino dominance), this is sufficient to 
wash out the order corrections. ^ The leading order results which we present here 
give the simple relations between the see-saw parameters necessary in order to obtain the 
phenomenologically successful LMA MNS matrix. Ultimately the masses and mixing angles 
in unified models must be calculated numerically, including the radiative corrections. The 
purpose of the analytic results we present here is to provide insight into the construction 
of unified models which must then be studied numerically. 

The paper is organised as follows. In section 2 we review the construction of the MNS 
matrix starting from general complex neutrino and charged lepton mass matrices, with 
equivalence of different parametrisations of the MNS matrix discussed in appendix We 
also discuss the charged lepton contributions to the MNS matrix, where the natural ex- 
pectation is that the charged lepton mixing angles are all small, and give an expansion of 
the MNS matrix to leading order in the small angles, with details relegated to appendix 
0. In section 3 we construct the successful neutrino mass matrix structures in a natural 
way using the see-saw mechanism with right-handed neutrino dominance. In section 3.1 we 



^By contrast in a bottom-up approach the order corrections may be considered [g^. However 
the effect of charged lepton mixing angles was not considered in jisj and this can significantly afi'ect the 
conclusions based on naturalness arguments. 
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derive results for single right-handed neutrino dominance with sequential sub-dominance, 
corresponding to a full neutrino mass hierarchy. In section 3.2 we derive results for single 
right-handed neutrino dominance with off-diagonal sub-dominance, corresponding to a par- 
tial neutrino mass hierarchy. In section 3.3 we derive results for off-diagonal right-handed 
neutrino dominance, corresponding to an inverted neutrino mass hierarchy, which we show 
is consistent with the LMA MSW solution once the effect of charged lepton mixing angles 
is taken into account. For each type of right-handed neutrino dominance we derive useful 
analytic expressions for neutrino masses and mixing angles in terms of the see-saw mass 
matrices, and show that for the hierarchical cases the LMA MSW solution gives a soft lower 
bound iC/esI ~ 0.1, just below the curent CHOOZ limit. The results obtained in section 
3 rely on analytic diagonalisation methods discussed in appendices ^ |^ In section 4 
we give two physical application of the results. In section 4.1 we consider Pf5Q^, and show 
that both hierchical and inverted hierarchical cases predict small f3f3oi, with \mee\ ~ 0.007 
eV within the sensitivity of future proposals such as GENIUS. In section 4.2 we discuss 
leptogenesis, and show that successful leptogenesis is possible if the dominant right-handed 
neutrino is the heaviest one, but the leptogenesis phase is unrelated to the MNS phases. 
Section 5 concludes the paper. 



2. Constructing the MNS matrix 

The charged lepton masses and the neutrino masses are given by the eigenvalues of the 
complex charged lepton mass matrix mfj^ and the complex symmetric neutrino Majorana 
matrix m^^, obtained by diagonalising these mass matrices. 



/me 

I (2.1) 

V 







/ mi 

m2 I (2.2) 



y m3 

where V^^ , V^^ , V^^ are unitary tranformations on the left-handed charged lepton fields 
El, right-handed charged lepton fields Er, and left-handed neutrino fields vl which put 
the mass matrices into diagonal form with real eigenvalues. 
The MNS matrix is then constructed by 

Umns = V''W^^ (2.3) 

In appendix |^ we discuss the relationship between different parametrisations of the MNS 
matrix. In this paper we find it convenient to parametrise the MNS matrix as 

Umns = U23U13U12 (2.4) 

which involves just three irremoveable physical phases 6ij. In this parametrisation the 
Dirac phase 5 which enters the CP odd part of neutrino oscillation probabilities is given 
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by 

S = Si3- 623 - 5i2. (2.5) 

where the phases are defined in appendix 

The MNS matrix is constructed in Eq. p.3| as a product of a unitary matrix from the 
charged lepton sector V^'^ and a unitary matrix from the neutrino sector V^^^ . Each of 



these unitary matrices may be parametrised by the parametrisation of V'^ in Eq. AJ. Thus 
we write 

Tjl^L rtfL Tjl'L TD^L Tyl^L Tjl^L /'o p,\ 



2 -"^23 -"'13-' 1 -"'12-' 3 
. tjEl pE 



pEL pE^ pEL pE^ pEL pE^ 
^2 -"'23 -"-13 ^1 -"-12 ^3 



(2.7) 



where the Euler angles and phases are defined as in Eqs.A.2-A.7 but now there are inde- 
pendent angles and phases for the left-handed neutrino and charged lepton sectors distin- 
guished by the superscripts and E^. As shown in appendix ^ the MNS matrix can 
be expanded in terms of neutrino and charged lepton mixing angles and phases to leading 
order in the charged lepton mixing angles which we argue must be small. By comparing 



Eq,B.12 to Eq.BJT we have: 



S23e 


-«<523 p 


" ^23^ 


^"23 — 


^23 '^23 e 


^°23 








0136 




"13 ^ 


^"13 — 


^13 '^23 ^ 




+ 






Sue 


-iSi2 p 


- ^12^ 


*''l2 -)- 


oEl „>^l p 
^23 ^12^ 




+ t 


()El fL fL JiQ -&f3 ) _ 

^13 ''12*23^ 


^12 ''23 ^-12^ 



^°12 



(2.8) 
(2.9) 
'(2.10) 



Clearly the large atmospheric and solar angles ^23 and 612 arise mainly from the neutrino 
sector with small corrections from the charged lepton angles. However ^13 receives impor- 
tant contributions not just from ^^3, but also from the charged lepton angles 9i2, and 
0^^ . In models where is extremely small, ^13 may originate almost entirely from the 

771 

charged lepton sector, since for example 6^2 niay be roughly equal to the Cabibbo angle 
in some models. Charged lepton contributions could also be important in models where 
^12 = very accurately, since corrections from the charged lepton mixing angles may 
allow consistency with the LMA MSW solution which requires 612 to be somewhat less than 
maximal. Such effects will be important for the inverted hierarchy model, for example, as 
we discuss in more detail later. 



3. See-saw models with right-handed neutrino dominance 

In this section we specialise to the case of the see-saw mechanism, with right-handed neu- 
trino dominance. A particular high energy theory will involve a charged lepton Yukawa 
matrix , a neutrino Yukawa matrix Y'^ and a right-handed neutrino Majorana matrix 
Mjiji. The Higgs vacuum expectation values (vevs) may be absorbed in the Yukawa matri- 
ces to give mass matrices mf^j^ = viY^ and m'^^ = V2Y'^, where we allow for two different 
vevs vi,V2 as is the case in supersymmetric models, while in the standard model vi = V2, 
and LR means that these are Dirac mass matrices in the left-right basis where the rows 
correspond to left-handed fields and the columns to right-handed fields. At high energies 
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these mass matrices rnfj^, m'^j^ and M/jjj are "born" in a particular basis which is defined 
by the particular unified theory. In general mf^^ and m'^^ are general complex matrices, 
while Mfiji is complex symmetric. 

In this "theory" basis the light effective Majorana neutrino mass matrix is given, up 
to an overall irrelevant sign, from the see-saw formula [16|,|29|, pO|]: 



m 



LL 



m 



LR^^^RR"''LR ' 



(3.1) 



ignoring the effects of radiative corrections. 

In order for our results to be maximally useful for top-down model building we must 
work in the "theory" basis defined by the high energy theory. In this basis all the mass 
matrices are in general off-diagonal, and in particular the charged lepton mass matrix will 
be off-diagonal, although according to our naturalness arguments we would expect it to 
yield small charged lepton mixing angles. Right-handed neutrino dominance should be 
applied in this theory basis, and questions of naturalness must be addressed within this 
defining basis. However it may happen that in specific theories certain simplifications 
naturally appear. For example the right-handed neutrino mass matrix may take one of the 



approximate forms in Eqs.1.3, 1.4, 1.5 



The general proceedure for diagonalising a general complex matrix is discussed in 
appendix |^. The strategy is to diagonalise the charged lepton and neutrino mass matrices 



mf^j^ and using, 



JD^L* U^L^ JD^L* t>El^ U^L^ JD^L* pEfi r/Efi r/Efi pEn jyEn pEfi 
^3 ^12 ^1 -"l3 ^21. ^2 ^LR^2 "^23 "^\3 ^\ ^12 ^3 



me 

rur 



pi/L* tjI^lT pUl^* tjI/lT ryi/i^T pi/i^* u pi^L* ni^L ni^L pi^L* r>i^L pi^L^ 
^3 -"'12 ^1 ^13 ^23 ^2 ''''LL^2 -"-23 -"-13 -'^1 -"-12 -'^3 



/ mi 



(3.2) 



m2 I (3.3) 

\^ TTls 

The specific case of a hierarchical neutrino mass matrix is discussed in appendix 0, while 
the inverted hierarchical case is discussed in appendix |E|. 

3.1 Full hierarchy from sequential sub-dominance 

In this sub-section we shall consider the case where the right-handed neutrino mass matrix 
takes an approximately diagonal form as in Eq jl.3| . If it is approximately diagonal, then it 
may easily be rotated to exactly diagonal form with real eigenvalues, to give 

/ X' \ 



diag 
RR 



0X0 

y yy 

In the theory basis we shall write the complex Dirac mass matrix as 

I a! a d 



(3.4) 



^LR 



b' b e 



(3.5) 
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Then using the see-saw formula |3 
symmetric neutrino mass matrix, 



(vahd for complex couplings) we find the complex 



' X' ^ X ^ Y 



\ 



a'b' 
X> 

X' ^ X ^ 



+ — + 

^ X ^ 



de 
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a'c' 
X' 
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ac 
X 


+ 


df \ 
Y 


e2 
Y 


b'c' 
X' 


+ 


be 
X 


+ 


ef 
Y 




c'2 
X' 


+ 


c2 
X 


+ 


fl 
Y / 



(3.6) 



The condition for single right-handed neutrino dominance with sequential sub-dominance 
was given as |24| 



Y XX' 
where x,y £ a,b,c and x',y' £ a',b',c', and now we are dealing with complex matrices 
we must consider the absolute value of the elements. For sequential subdominance we can 
essentially ignore the contributions from the right-handed neutrino with mass X' , so that 
the problem reduces to the simpler case of two right-handed neutrinos [p3] with 



m 



LL 



X 
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^ Y 


ab 
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de 
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ac 
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+ 
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^12 






62 
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e2 
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be 
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Y 










^23 










c2 
X 


+ 


^) 




K 


^13 


"^23 


"^33 



(3.8) 



where we have made contact with the notation of Eq.D.l. In order that Eq.3.8 resembles 



the leading order form of the matrix in Eq.|Ll] we require 

\d\ « \e\ « I/I 



(3.9) 



as first shown in |22|, although in the complex case here we have taken the absolute values 
of the complex couplings. It is straighforward to diagonalise the matrix in Eq. |3.8| by 
following the proceedure in appendix According to Eq. D.l there is a complex phase 



associated with each element of the mass matrix. In our previous results the effect of 



phases was ignored |24|. Here we include them, and discuss their effect on our previous 
results. 



We first focus on the lower block of Eq.3.8 



771,22 



23 

"^23 "^33 



^ _L e! bc^£f 

X ^ Y X ^ Y 

^ _L £/ i£ I /! 

X ^ Y X ^ Y 



(3.10) 



This matrix may then be diagonalised by the ^ re-phasing followed by the 23 rotation 



m 







m'o 



22 



R 



23 



777^3e-*(<^: 



2j(/>2^ 



2 



777236 



m 



^23 



(3.11) 



33": 



The 23 mixing angle may readily be obtained from Eq. p.3| , which is accurate to order ^13. 
We shall find that in the hierarchical cases ^13 ^ 7772/7773 and that this bound must be 
saturated for the LMA MSW solution. Therefore in this case the results will be accurate 



to order 7772/7773. From Eq. D.3 we find 



tan 26*23 



Y 



-2ich% e2 



) 



(3.12) 
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where 



be 

a — X 



^ X 



X 



In the leading approximation, in which the corrections of order are neglected, the phases 
arise solely from the complex Dirac couplings e = |e|e*'^'=, / = \f\e^'^f, and we have the 
elegant leading order result 

|e| 



tan 



(3.14) 



where the phases are fixed by 



6- 



(3.15) 



This leading order result was first written down without phases in |2^. The result in 

^23 



Eq,3.14 demonstrates that the phases on e, / are not important in determining 



leading order, only their absolute values matter. 



It is instructive to take the determinant of both sides of Eq.3.11 



77122^3 



771226 



33 



-2i(j>„ 



(m2^3)2e-2^(*2- 



-2i{0g+<^^^) 



ce] 



XY 



(3.16) 



where the leading order terms proportional to 1/y^ have cancelled in constructing the 
determinant. This is no accident: it happens because the determinant vanishes in the limit 
of a single right-handed neutrino. The fact that the determinant is small implies that 
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(3.17) 



The natural origin of Eq,3.17 is crucial both in obtaining a natural hierarchy m2 <C 
and in obtaining a large solar angle. The naturally small 23 sub determinant is the main 
consequence of single right-handed neutrino dominance, as emphasised in |23]. 



If we take the trace of both sides of Eq. 3.11 , using Eq. 3.17 , then we may obtain the 
third neutrino mass to leading order 



m 



/ i(20--202^) 



Y 



(3.18) 



7773 



Y 



(3.19) 



A more accurate expression for 7773 including the 7772/7773 corrections can be readily obtained 
from Eq. p.^ , using Eqs. 3.10| , 3.12| . However the result is not very illuminating, and so is not 
worth displaying explicitly, although it may be readily constructed if required. The same 
comments apply to the remaining parameters, which we only display explicitly to leading 
order. 



The leading order result for 777,22, then follows from Eqs. 3.18 , 3.16 , 3.14 , 3.15 



"1-22 



-2i</.2 



-23' 



■^23 



-rf)\2 



X 



(3.20) 
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The leading order results for m^g; "^13 follow from Eq. D.9 , which we write here as 



m 



m 



12 

u 
13 



''23 



(3.21) 



'-13': 



We find to leading order using Eqs j3.21| , p^j3.14| , ^.15 



X 



(3.22) 



m 



13 



(3.23) 



Note that the 1/Y terms have cancelled to leading order in Eg. 3.22 as in the real case 



The 13 neutrino angle is given from Eq.D.lC 



h3 



m' 



13 



(3.24) 



where the leading order form for is given in Eq.3.23, and for m'^ in Eq,3.18. In the 
large d limit 

\de\ \df\ \ab\ \ac\ 

— »^,^ (3.25) 



we find the simple leading order result |24], 



h3 



\d\ 



VM' + \f\' 



with the phase (/ig^ fixed by 



(f>2^ ^ (l)e- ff^d 



Ignoring phases, this was the result previously quoted 
d limit 

\de\ \df\ \ab\ \ac\ 



Y ' Y 



X ' X 



(3.26) 
(3.27) 

On the other hand in the small 

(3.28) 



we obtain the leading order result 

1 



,::.a(4f6 + c^tce^(<^^^-'^/)) 



Qt'L ^ " ^-i't'o "-V-^23'^ "I" ^-23' 

^ Y 



(3.29) 



where the phase i;^2^ is fixed by Eq.D.ll 



After the 13 rotation the neutrino mass matrix is in block diagonal form 





^11 


m^2 







™12 


"i22 





v 








m 



(3.30) 



where the leading order form for rh\2 is given in Eq.3.22, for 77122 ™ Eq.3.20. We need to 
find from Eq.D.lS 



rrin 



(3.31) 
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Using the leading order form for m^'f in Eg. 3.23 , and for m'^ in Eg, 3. 18 , and m^^ in Eg.| 
we find 



m 



11 



a 

'x 



2d 



X 



vw+w 

where the leading 1/Y dependence in m^^ has cancelled, as in the real case [24]. 



(3.32) 



In the leading order, valid up to corrections of order 7712/7713, we may ignore the second 
term in Eg. 3.32 since this cannot be larger than order , as shown in Eg.3.2f:. Then the 



upper block which remains to be diagonalised in Eg.3.3C is given from Egs.|3|, |]2|, |330| 
and after the phase transformations P^^ it is: 



AB 
AB B^ 



where 



A 



B 



The leading order form in Eg. 3.33 clearly has a vanishing determinant 

0. 



777-'^ 7712 



(3.33) 



(3.34) 



(3.35) 



This is no surprise since seguential sub-dominance corresponds approximately to the case 
of two right-handed neutrinos, and hence the lightest neutrino mass is approximately zero, 
as in the case of a full hierarchy. The situation is analagous to the case of the lower block 
with single right-handed neutrino dominance discussed above. The second neutrino mass 
is simply given at leading order from the trace of Eg. 3.33| , 



m'2 



A^ + B'^ 



The leading order 12 neutrino mixing angle is given from Eg. p.l4| , |3.35 



ae 



(3.36) 



(3.37) 



B (c^t?,-4tce'(*^-'^?)) 
which is analagous to the real result in [^]. The phase x'^^ is fixed to give a real 12 angle 

|6| sin((/>'f,) « Sat |c| sin((^^) 



as in Eg.D.15 

and then 



tan 



C2^\h\ COs(0'f,) - 83^ |c| cos((/)^) 



where 



K = 0fe - 0a - 02 
b'^ = (l)c - (t>a + 4>e 



X 



(3.38) 
(3.39) 



(3.40) 
(3.41) 



Eg, 3.39 shows that a large solar angle 0^2. ~ '^/^ reguires |a|/(|6| cos(0^) — |c| <io?>{j^'^') ~ 1 
which is the complex analogue of the real condition \f2al{}j — c) ~ 1 derived in 
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The effect of phases is now quite important since the denominator of Eq. 3.3£ may involve 
phase-dependent cancellations. In the absence of such cancellations, the basic physical 
requirement for a large mixing angle, namely that the sub-dominant Dirac coupling a be 
of order b,c remains as in |24|. 

From Eqs. |3l8|j33^j335 , we see that sequential sub-dominance has generated a full 
neutrino mass hierarchy 

mi <C m2 <C ms (3.42) 



and hence 



Am|2 



m 



35 



Am2i 



rrio- 



(3.43) 



Sequential sub-dominance naturally leads to the LMA MSW solution by assuming a mild 



hierarchy of couplings in Eq.3.7. 

We now discuss an order of magnitude lower bound on the 13 neutrino mixing angle. 
Clearly such a bound can only be order of magnitude since it is always possible to arrange 
vanishingly small values of this angle by tuning parameters. However, in the absence of 
tuning, it is possible to make some general statements about the expected magnitude of 
this angle. The question is therefore, in the absence of unnatural cancellations, how small 
can the 13 neutrino angle be? Clearly from Eqs. 3.24, 3.23] , in the absence of tuning, the 



smallest values of correspond to setting d = 0, and hence the limit will be saturated 



by Eq.3.29. The 13 neutrino angle in Eq.3.29 cannot be made arbitrarily small by setting 
a = since we have just seen that a large solar angle requires a be of order b, c. Therefore 



the numerator of Eq j3.29| is expected to be of a similar order of magnitude to m2 in Eq j3.36 
and for a large solar angle we therefore obtain the order of magnitude lower bound 



h3 I 



|m2/m3| 



(3.44) 



From Eq. B.12, since the charged lepton angles are all required to be small due to the 



naturalness arguments , (^13 
bound on the MNS element, ^ 



013, and using Eq. p.43 , we may elevate Eq. 3.44 to a 



\Ue3\^> \ Ami J Am] 



2 I 

■32 1 



(3.45) 



The bound is roughly proportional to tan 0^2) by eliminating a. Interestingly Eq.3.45 
predicts that for the LMA solution \Ue3\ should be just below the current CHOOZ current 
limit, which raises the prospect that it could be measured at the forthcoming long baseline 
(LBL) experiments. We emphasise again, however, that the bound is a soft one since it is 
based in the premise of there being no cancellations in the construction of the MNS matrix, 
and so is really only approximate to within an order of magnitude. 



3.2 Partial hierarchy from off-diagonal sub-dominance 

We now turn the case where at high energies the upper block of the right-handed neutrino 
mass matrix has the approximate off-diagonal form as in Eq.lL^. By small angle rotations 



^Similar bounds have been obtained previously in somewhat diflerent frameworks in ||3^,p^ 
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it may then be rotated into the form below with real mass parameters X, Y, 

f X 



M 



off— diag 
RR 



X 

y 



(3.46) 



Using the same notation for the Dirac couplings as in Eq.3.5, the see-saw formula in Eq.3.1 
then gives in this case, 



/ 2aa' I d?_ a'b , aV , de a'c i ac' • df \ 



m 



LL 



Y X ^ X ^ Y X ^ X 

2bb' , ^ I bc^ I 



X 



\ 



Y 

_ _ 
X ^ X ^ Y 

X ^ Y / 



(3.47) 



The condition for single right-handed neutrino dominance with off-diagonal sub-dominance 

\e\\f\\ef\ 



was given as |24| 



\xx 



Y » ^ (3-48) 

where x E a,b,c and x' G a' ,b' ,c' , and now we are dealing with complex matrices we must 
consider the absolute value of the elements. 

we 



As before, for the neutrino mass matrix to resemble the leading order form in Eq.1.1 



also require the condition in Eq.3.9. Then the discussion of the leading order 23 neutrino 
mixing angle follows exactly as in the sequential sub-dominance case, with the lower block 
replaced by 



^22 "^23 



'-23 



2bb' 
X 



b'c 



+ 



be' , ef 



X 



b'c I be' I ef 
X ^ X ^ Y 



X 
2ec' I 
X ^ Y 



Y 



From Eq.D.3 we find the analagous result to Eq. f3.12| , 



tan 26*23 



where now 



ye ,bci 
X ^ X 

£/ 
Y 



4 



-2i<t>'iP. 
Y 



^-2ir,^) 



„-2i(h'< 2ee' 

e -^i — 



-2i4>Y^ 2hy_ 



Y 



e Y 



(3.49) 



(3.50) 



(3.51) 



In the leading approximation, in which the corrections of order are neglected, we obtain 
the same results as in Eqs.3.14,3.15, 



tan 9<^Q ~ 7T7- 
I/I 

The determinant of Eq.3.11 with Eq. 3.49| gives 

-2»(</,!;+<^;;^) (bf - ce){b'f - c'e) 
XY 



2e' 



(3.52) 



(3.53) 



where again the leading order terms in the determinant have cancelled due to single right- 
handed neutrino dominance. The trace of Eq.3.11 with Eq.3.49 gives the same leading 
order result for mg is as in Eq. 3.1^ , 

P + l/P] 



„i(2</,--2</,2^)! 



Y 



(3.54) 
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Hence from Eqs.3.55 , 3.54, 



m 



22 



2e 



~2j0: 



X 



(3.55) 



As before, henceforth we display only the leading order results explicitly, with the correc- 
tions 1112 /m^ readily obtainable from our formalism if required. The leading order results 



for 17112, "^13 follow from Eq.3.21. We find to leading order using Eqs,3.21,3.47 



X 



(3.56) 



m 



13 



+ 



X 



+ e 



i(j>, 



d^/W+W 



Y 



(3.57) 



Again the leading order 1/Y terms have cancelled in Eq.3.56 



The 13 neutrino angle is given from Eq.D.lO where the leading order form for m^g is 
given in Eq. |3.57 , and for rrig in Eq.3.54. In the large d limit we find the same result as in 



Eq,3.26 



h3 



\d\ 



vWTW 

On the other hand in the small d limit we find the leading order result 



^3 ^ m' 



X 



+ 



-23' 

X 



(3.58) 



(3.59) 



where the phase (j)2^ fixed by Eq.D.ll 



After the 13 rotation the neutrino mass matrix is in block diagonal form as in Eq,3.30 
Using the leading order form for m^k^ in Eq.3.57, and for m'^ in Eq.3.54, and m^^" in Eq,3.47 
we find 



11 



2aa' 



2d 



V\e\' + \f\' 



(.^3-6 + c^fce^(^^-'^P) , a(4f6' + 4fc'e^ 



X 



+ 



X 



(3.60) 

where the leading 1/Y dependence in m^^" has cancelled, as before. After the phase trans- 
formations the upper block which remains to be diagonalised is, approximately. 



2AA' A'B + AB' 
A'B + AB' 2BB' 



where 



A' 



B' 



X 



(3.61) 



(3.62) 



and A,B are as in Eq. 3.34 . The determinant of Eq. 3.61 is now non-zero due to the two 
sub-dominant right-handed neutrinos. 



m 



^ -{A'B-AB') 



l\2 



(3.63) 
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The trace of Eq.3.61 is, 



m[ + 2{AA' + BB') 



The leading order 12 neutrino mixing angle is given from Eq.D.14,3.61 



tan 



tan (eir'' + ^12^) 



where the angle 9^2 is given by 



A 
B 



ce 



and is given by 



tan 9'^2 



B' 



(3.64) 



(3.65) 



(3.66) 



(3.67) 



The phase x'^^ is fixed by a similar proceedure to that described previously. 

From Eqsf]5|Jp|,p|,|l| we see that off-diagonal sub-dominance has generated a 
partial neutrino mass hierarchy 

(3.68) 



vn\ ^ 7712 ^ ms 



and hence 



m 



3' 



Am: 



21 



2 



(3.69) 



Off-diagonal sub-dominance also naturally leads to the LMA MSW solution by assuming 
a mild hierarchy of couplings in Eq. 3.4^ . 

Analagous to the sequential sub-dominance case, there is a lower bound on 9^^ arising 
from small d limit limit Eq. 3.59| . If we assume no accidental cancellations then the square 



of the numerator in Eq. |3.59| is of the same order as the product of eigenvalues in Eq j3.63 



Since the lightest two eigenvalues are of similar magnitude in this case, we may deduce a 
lower bound analagous to Eqjslsj 



(3.70) 



3.2.1 The pseudo-Dirac limit 

A limiting case of off-diagonal sub-dominance is when the lightest two neutrinos form 
an approximate pseudo-Dirac pair, which happens when the terms in the upper block in 



Eq,3.61 satisfy either 



\J^V\B\ > \A\AB'\ 



or 



\A\,\B'\ > l^'U^I 
The conditions for this happening are, respectively, using Eqs.|]6|,|^, 



tan0j'fi<Cl, tan0^f>l. 



or 



tan 



'12 



> 1, i?m9\}i<^\ 



(3.71) 
(3.72) 

(3.73) 
(3.74) 
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In either of the two Hmiting cases, from Eq. ^.65 , we will have an almost maximal neutrino 
contribution to the solar mixing angle, 



'12 



7r/4 



(3.75) 



and, since the trace of the upper block is very small, opposite sign mass eigenvalues. 



-rrio 



(3.76) 



both of which are characteristic features of a pseudo-Dirac neutrino pair. In the pseudo- 
Dirac case we clearly have 



Am2i ^ 1712 



(3.77) 



The hierarchy in Eg. 3. 77 implies that the pseudo-Dirac case is applicable to the LOW 
solution rather then the LMA MSW solution, so is not of interest to us here. 

Finally, the 13 neutrino mixing angle in the pseudo-Dirac case is given in the small d 
limit from Eq,3.59, using Eqs. 3.71, ^.72 , 



^13 



m. 



m 



(3.78) 



Hence in the pseudo-Dirac case we have the bound analagous to Eq. |3.44 

l^ifl > Ims/msl (3.79) 
This bound differs from that quoted in [p^], since it is not singular for 612 = '/r/4. 



3.3 Inverted hierarchy from off-diagonal dominance 

We now discuss the case of off-diagonal right-handed neutrino dominance which gives 
rise to a neutrino mass matrix as in Eq jl.2| corresponding to an inverted neutrino mass 
hierarchy with opposite sign mass eigenvalues for the heavier pseudo-Dirac neutrino pair. 
The mechanism is based on the off-diagonal form of the right-handed neutrino mass matrix 



as in Eq.3.46, which after the see-saw mechanism gives rise to the physical neutrino mass 



matrix in Eq.3.47, 



LL 



/ 2aa' I cP a'b , aV i de a'c i ac' , df_ 

ef 



X 



Y X ^ X 

2bb _ 
X ~^ Y 



Y X ^ X 
I ef. ^ j_ be' 



V 



2cc' 
X 



^ Y 



(3.80) 



However, instead of applying the condition for single right-handed neutrino dominance as 



in Eq.3.48, we shall now require the pseudo-Dirac pair of right-handed neutrinos with mass 
X to dominate in such a way as to give the leading order form of the neutrino mass matrix 



in Eq. |l.2| . The two alternative conditions for this are either: 

\a'b\ 



X 



\a'c\ 

IT 



> 



\a'a\ 



\x'x\ 



X ' X 



yy 

Y 



(3.81) 
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where x' G 6', c', and x £ a,b, c, and y, y' £ d, e, f. 
Or: 



\ab'\ 



lac'l laa'l Ixx'l \vv'\ 



(3.82) 



X X " X ' X ' Y 

where x E b,c, and x' S a',b',c', and y,y' G d,e,f. In the following we shall consider the 
first alternative in Eq.3.81 for definiteness. 

The proceedure for diagonalising the form of the mass matrix in Eq. 3.8C| with the 
inverted hierarchy conditions in Eq.3.81 is outlined in appendix ^. We first perform the 



re-phasing as in Eq.p.2|. Then we determine the 23 neutrino mixing angle OXi from Eq.E.5 



tan 



'23 



(3.83) 



where 



ac' I df 
X ^ Y 

a'c 

X 



ab' I de 
X 



The leading order result, neglecting the corrections ii is then 



(3.84) 



tan^^f ~ 



(3.85) 



where the relative phases are fixed by 



(3.86) 



We find at leading order 



X 



(3.87) 



Henceforth we only give the leading order results explicitly, as in the previous cases. The 
lower block elements are given by 



m 



23 ^23 ^23 ^ 



23 



X 



Y 



{31 



m. 



22 



X 



+ 



Y 



"L I Sno e — c: 



-j2</>2^ l'^23' 



Y 



(3.89) 
(3.90) 



We next perform the small angle 13 rotation in Eq.E^. From Eqs.E.ll,3.87,3.8S we deter- 



mine 



'13 ' 



H3 



m. 



23 



m 



(3.91) 



12 



The requirement that is real fixes the absolute value of the phases 



'2 1 v^3 



-18- 



After the P^^ re-phasing the upper 2x2 block which remains to be diagonahsed is of 
the form 



m 



(3.92) 



where from Eq.3.87,3.8£, 



m 



12 



1 ,2aa' cP 
7(^ + 



m 



X 



22 



m 



(3.93) 
(3.94) 

(3.95) 



12 



Eq.3.92 is of the pseudo-Dirac form. The trace and determinant of Eq.3.92 reveal the 
two heavier eigenvalues m'^ , m'2 , 



m'l + m'2 = {5 + e)m' 
m'im'2 = (5e — l)(m')^ 



Thus we have 



-m'l 



m'2 



m 



Am 



32 



\m'2\' 



Amli w 2\m'2\'^Re{6 + e) 
We find the neutrino contribution to the 12 mixing angle by diagonalising Eq. 3.92| , 

tan^r2 



(3.96) 
(3.97) 

(3.98) 
(3.99) 
(3.100) 

(3.101) 



which is almost maximal due to the smallness of 5, e or equivalently the pseudo-Dirac form. 
From Eqs. |3.10C|j3^ we find 



Aml2\ 



2Re{5 + e) 



From Eq.3.101 we find 



tan (9if 



(3.102) 



(3.103) 



These relations together with Eq.3.91 give important constraints for the inverted hier- 
archical case. By comparing Eq. 3.95 to Eq. 3.9l| it is apparent that, in the absence of 
cancellations, O^^ ~ e. The parameter 6 is completely independent of O^k, and e. Thus for 
example we can have O^i^ = e = 0, while at the same time having 5 7^ corresponding 
to -R 7^ 0. Thus, unlike the hierarchical case, there is no lower bound on O^^. On the 
other hand there is an upper bound on , in the absence of cancellation effects, which is 

and 5 we find 



saturated when 5 is not large compared to e ~ . For all choices of e 

\el^\ < ll-tan^^'fl ~i? 



(3.104) 
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Since for the LMA MSW solution we know that R ~ 10"^, Eq j3.104| imphes that the 
neutrino mixing angles = 7i"/4 and 9\i^ = to high accuracy. However the LMA MSW 
solution is not consistent with maximal mixing and has a preferred mixing angle 612 ~ 7r/6, 



and for this reason there have been claims in the literature (see for example [281) that 
the inverted hierarchy case is inconsistent with the LMA MSW solution. However this 
conclusion ignores the contribution to MNS angles arising from the charged lepton sector, 
as was originally pointed out in |26|. It is clear from Eqs.p.9, 2.10| that even if 0'^2 — '^/^ 
and = that a reasonable charged lepton contribution, for example a Cabibbo-like 
contribution to ~ 0-2, can result in 612 in the LMA MSW range and 6*13 close to the 
current CHOOZ limit. 

4. Physical applications 

4.1 Neutrinoless double beta decay 

There has been a recent claim of a signal in neutrinoless double beta decay correponding 
to |mee| = 0.11 — 0.56 eV at 95% C.L. However this claim has been criticised by two groups 



1 33], 1 34] and in turn this criticism has been refuted ]35]. Since the Heidelberg-Moscow 
experiment has almost reached its full sensitivity, we may have to wait for a next generation 
experiment such as GENIUS to resolve this question. 

From the theoretical point of view, in the natural hierarchical and inverted hierarchical 
models favoured here, corresponding to the leading order mass matrices in Eqs jl.l] , |1.2| , it 
is immediately clear that very small values of \mee\ *C 0.05 eV are expected, since both 
matrices have leading order zeroes in the 11 position. Since we argue that these are the only 
two natural leading order forms, we can immediately predict that the signal for neutrinoless 
double beta decay should be below the sensitivity of the Heidelberg-Moscow experiment. 
The next question is whether a more sensitive experiment such as GENIUS is likely to see 
a signal? 



This question has been fully analysed in |33], using the formula 



E(f^er')Hl (4.1) 



\mee 

3 

and the general conclusions are that for the hierarchical models ]mee] = 0.0005 — 0.005 eV 
while for the inverted hierarchical models \mee\ = 0.01 — 0.057 eV [^3|. It is interesting to 
re-examine these expectations from the point of view of the equivalent formulation 



me 



\iV^'ml^V^^')n\ (4.2) 



If we neglect the charged lepton contributions then \mee\ is just given by ](?tt.^^)ii] and 
we can read-off the prediction for this in the different models. Of course the contributions 
from the charged lepton sector can be very important, as we have emphasised, however 



what can look like a cancellation conspiracy in the formulation in Eq. [4.l| can be readily 



understood in the formulation in Eq.4^, so the alternative formulation can also be quite 
instructive, and we use it here as a simple application of our analytic results. There are 
three case to consider and we discuss each in turn. 



-20- 



We shall show that in all the cases considered, full hierarchy, partial hierarchy, and 
inverted hierarchy, where the hierarchies result naturally from right-handed neutrino domi- 
nance, that neutrinoless double beta decay should have a value of about of \mee\ ~ 0.007 eV 
assuming the LMA MSW solution. Such a value is below the sensitivity of the Heidelberg- 
Moscow experiment, but is within the range of future proposals such as the GENIUS 
experiment. 

(i) Full hierarchy from sequential subdominance 



In this case, ignoring the charged lepton rotations, from Eq.3.8 we find 



l^eel f« |(mii)ii| w 1^ + ^1 (4.3) 



Ignoring the contribution from y- which is always small even in the large d limit in Eq,3.25 



and using Eq.S.SC we find 



mee\ ~ m2 « V Am|i ^ 0.007eV (4.4) 



which slightly exceeds the upper end of the quoted hierarchical range |33] 
(ii) Partial hierarchy from off-diagonal subdominance 



Again, ignoring the charged lepton rotations, from Eq. 3.47 we have 



|mee| ~ |(mii)ii| PS 1^ + yl (4.5) 



Again, ignoring the contribution from ^ which is always small, and using Eq.3.64 we again 
find 



\mee\ ~ ms PS y Am|i ^ 0.007eV (4.6) 

In the pseudo-Dirac limit the term would be very small, so much smaller values of 
|mee| would be expected. However, as we discussed earlier, the pseudo-Dirac limit is not 
relevant for the LMA MSW solution. 



(iii) Inverted hierarchy from off-diagonal dominance 



Ignoring the charged lepton rotations as before we find, from Eq.3.80, 



Im^el « |(m^i)n| « |^ + — I (4.7) 



Using the results of the previous section with Eq.4.7 we find 



\5.m2\ ^ \Rm2\ ^ 0.01m2 < 0.01a/ Am^a ^ 0.0005eV (4 
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The result in Eq. |4.8| neglects the effect of charged lepton mixing angles. However we saw in 
the previous section that the inverted hierarchy case in inconsistent with the LMA MSW 
solution unless the effect of charged leptons is considered. Therefore one would expect the 
result in Eq.^I^ to be considerably affected by charged lepton mixing angles in this case. 
As an example let us consider the effect of the 12 charged lepton mixing angle on the value 



of \mee\- From Eq. 4.2 we find a contribution from 6^2 ■, 



\mee\ ^ |^|A/|Am|2| ~ |-P-|0.008eV (4.9) 
V2 V tic 

where we have scaled the charged lepton angle by the Cabibbo angle Oc- Thus, if the 
inverted hierarchy case is to achieve consistency with the LMA MSW solution, we would 



expect a contribution to neutrinoless double beta decay to be given by Eq. 4.£. 



4.2 Leptogenesis 

Leptogenesis pH], [37] provides another example where a correct treatment of phases is 



crucial to the physics, and this provides another example of the application of our general 
treatment including phases considered in this paper. An analytic treatment of leptogenesis 
in the the framework of single right-handed neutrino dominance was considered earlier p8| . 
In this subsection we highlight how the more general treatment of phases considered here 
impacts on our previous analytic estimates. The numerical results previously presented 



1 38] are of course completely unchanged. 

According to leptogenesis the baryon asymmetry of the universe is given by 

Yb = (4.10) 

a — 1 g* 

where a is a number of order unity arising from sphaleron effects which convert the lepton 
number to baryon number, d is a dilution factor which is crucial to the success of leptoge- 
nesis since it can be many orders of magnitude below unity if the right-handed neutrinos 
are not produced efficiently enough or if they do not satisfy the decay-out-of-equilibrium 
condition strongly enough, and g* ~ 10^ represents the number of effective degrees of free- 
dom. In the analytic discussion here we shall confine our attention to ei which is the CP 
decay asymmetry parameter of the lightest right-handed neutrino, and is the basic seed of 
lepton number violation in the leptogenesis scenario. 

^ r(iVm ^ + H2) - r(< ^ L/ + Hi) 
T{Nri ^ L,- + H2) + T{n\^^ ^ L/ + ^t) 

Assuming a hierarchy of right-handed neutrino masses Mi ^ M2 <^ M3, one has, approx- 
imately. 



mnMn)^]') (^) (4.12) 



Leptogenesis therefore depends on following Yukawa combinations, or, equivalently, Dirac 
mass combinations 

{iyi:R)^Yi:j,)^, = J.((m^^)t^^^)i, (4.13) 
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The combinations in Eg, 4. 13 have the property that they are invariant under a change 
of charged lepton basis, so they may be calculated either in the original theory basis or 
in the diagonal charged lepton mass basis. However the combinations are sensitive to a 
basis change of the right-handed neutrinos. For the present discussion we shall assume 



that in the theory basis the right-handed neutrinos are diagonal as in Eq.1.3. We shall also 
assume right-handed neutrino dominance with sequential subdominance. Within this class 
of model, leptogenesis provides a means of discriminanting between different subclasses 
of sequential subdominance, namely between the case where the dominant right-handed 
neutrino is the heaviest {Y = M3) and the case where it is the lightest {Y = Mi), where 
we assume Mi <^ M2 <C in both cases. We found that leptogenesis actually prefers 
the case Y = M3 where the dominant right-handed neutrino is the heaviest The 
reason is that leptogenesis is only marginally consistent with the gravitino constraint on 
the reheat temperature Tr < 10^ GeV, and the values of ei for the case where the dominant 
right-handed neutrino is the lighest one are suppressed by a factor of 7712/7713 relative to 
the case where it is the heaviest, and this suppression is enough to make it quantitatively 



unnacceptable [38|. We now consider the two possibilities explicitly. 



(a) Sequential subdominance with Y <ti X <^ X' 



In this case we need to re-order the matrices involving right-handed neutrinos so that 



the first column corresponds to Y, and so on. Then the combinations in Eq. [4.13| can be 
written as 

^ < Zl\zi > < Z\\Z2> < Zl I Z3 > \ 
((777^^)^777^^) = < Z2\zx > < Z2\z2 > < 2:2^3 > (4.14) 
\< Z3\zi > < Z3\Z2 > < Z3\Z3> J 

where we define three complex vectors \zi >= {d,e,f), \z2 >= {a,b,c), \z3 >= {a',b',c'). 



Clearly the combinations in Eq.4.13 only depend on the first row of the matrix in Eq,4.14 
and hence only on two phases, namely that of < zi\z2 > and < zi\z3 >. ^ From the 
sequential subdominance conditions in Eqs. |3.7| , we find that one term dominates 



3 Y lra{< zx\z2>^) 



X^'Kvl X 



< Z\\Z\ > 



(4.15) 



so that leptogenesis only depends on one phase combination. Using Eq.3.8, and the leading 
order results this may be expressed as 



Y 



Wttv^ X 



Im ( e 



b + C23 ce 



Using Eqs. p^^ we find the order of magnitude result 



~ 2 X 10"^— sin(20i2' 
7773 



^vlOiOGev) 



(4.16) 



(4.17) 



*In general a third phase wiU enter leptogenesis via < 22 1^3 >, but in the approximation that lepton 
number arises from the decay of the lightest right-handed neutrino it is not relevant. 
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where 612 = arg < zi\z2 >, which quantitatively turns out to be too small due to the ^ 
suppression factor |38|. 

(b) Sequential subdominance with X' <^ X <^ Y 

The case where the dominant right-handed neutrino is the heaviest requires no re- 
ordering of the mass matrices, and we have Eq. |4.14| with \zi >= (a', b' , c'), \z2 >= (a, b, c), 
1^3 >= {d,e,f) (if Ml = X',M2 = X,M3 = Y). From the sequential subdominance 
conditions in Eqs.|3.7|, we find the dominant term is now 



_ 3 X' Im {< zi\zs y"^) 
^ IGttv^ Y < zx\z\ > 



(4.18) 



so that leptogenesis now depends on the < ^il-zs > phase combination. Using Eq.^^, and 
the leading order results this may be expressed as 



X' |e|2 + |/|^ 



\^-Kvl Y la'P + |6'|2 + |c' 



^2 

Using Eqs.p^, we find the order of magnitude result 



-^Im (e2^<^^^(4t(6')* + c^t(c')*e^^*^"^^¥) (4-19) 



) ~ 2 X lO-Sin(20i3) (^^T^) (4.20) 

where <\)\-i = arg < zi\z3 >, which quantitatively turns out to be just about acceptable 
since it does not suffer from the ^ suppression factor of Eq. 4.17 |]38|| . 



Note that the leptogenesis phase (f>i3 in this case involves the phases of the Dirac 
couplings to the lightest right-handed neutrino of mass X' . These couplings and phases 
are completely irrelevant for constructing the MNS matrix. 

5. Summary and Conclusion 

The latest solar and atmospheric neutrino data strongly point towards a minimal interpre- 
tation consisting of three active neutrinos whose mixings are described by a LMA MNS 
matrix, corresponding to the mixing angles 623 ~ 7i"/4, 9i2 ~ vr/6, ^13 ^ 0.2. We have 
argued on naturalness grounds in favour of a neutrino mass matrix whose leading order 
form is either Eq.^]^, corresponding to a neutrino mass hierarchy, or Eq.1.2 corresponding 



to an inverse mass hierarchy with the heavier neutrinos comprising an approximate pseudo- 
Dirac pair. We have constructed the LMA MNS matrix to leading order in ^13 including 
the neutrino and charged lepton mixing angles and phases, the latter playing a crucial role 
in allowing the inverted hierarchy solution to be consistent with the LMA MSW solution. 



We then showed how the neutrino mass matrices with leading order forms in Eqs.1.1,1.2 



may be constructed naturally from the see-saw mechanism with right-handed neutrino 



dominance. For the hierarchical form of the mass matrix in Eq.Ll, single right-handed 
neutrino dominance ensures that we obtain a natural neutrino mass hierarchy m2 ^ 
with no tuning, not even at the level of 10% accidental cancellations, which would otherwise 
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be required. Although such a tuning is quite mild, it would mean that the low energy results 
would be subject to rather large radiative corrections, and therefore that the low energy 
neutrino masses and mixing angles are not directly related to the high energy ones. By 
contrast in technically natural theories based on single right-handed neutrino dominance, 
the radiative corrections are expected to be quite small, at the level of a few per cent p5|. 



For the inverted hierarchy case, where the leading order form of the mass matrix in Eq.1.2 



arises from off-diagonal right-handed neutrino dominance, the radiative corrections are also 
smah [|6|. 

With right-handed neutrino dominance because the radiative corrections are small the 
low energy observables provide a direct window into the high energy parameters of the 
theory. In such theories it is therefore meaningful to relate the high energy see-saw param- 
eters to the low energy neutrino observables by simple analytic relations which ignore the 
radiative corrections. We have done this for three types of right-handed neutrino domi- 
nance corresponding to a full hierarchy from single right-handed neutrino dominance with 
sequential sub-dominance, a partial hierarchy from single right-handed neutrino dominance 
with off-diagonal subdominance, and an inverted mass hierarchy from off-diagonal domi- 
nance. In each case we have derived analytic expressions for neutrino masses, mixing angles 
and phases in terms of the parameters of the see-saw matrices. Our results are accurate to 
leading order in ^13, which is sufficiently accurate bearing in mind that such expressions 
ignore the effects of radiative corrections. The leading order expressions are particularly 
instructive in guiding the construction of see-saw models, and these should be particularly 
useful for unified model building. This analysis extends the range of the analytic results in 
|24],||2^ to include the important effects of phases and charged lepton mixing angles which 
are often overlooked in the literature. 

It is worth highlighting how the large atmospheric and solar mixing angles arise for 
the different types of right-handed neutrino dominance. In all cases we write the neutrino 
complex Dirac mass matrix in the common notation of Eq.3.5 where the columns of this 
matrix represent the couplings to the different right-handed neutrinos to the right-handed 
neutrinos with Majorana mass matrices given in Eqs jl.3|jL5| . Without loss of generality the 
columns of the matrices may be permuted corresponding to a re-ordering of right-handed 
neutrinos. 

For case of a full neutrino mass hierarchy mi <C ?n.2 *S with the right-handed 
neutrinos having the texture in Eq.1.3 with the right-handed neutrino of mass Y giving the 
dominant contribution and the right-handed neutrino of mass X giving the sub-dominant 
contribution, and the right-handed neutrino of mass X' giving negligible contributions, 
the large atmospheric angle then arises from Eq. 3.14| , and the large solar angle arises from 
Eq j3.37 . The fact that the solar angle is given, up to phases, by the ratio of the subdominant 
coupling a to the difference of subdominant couplings b — c, was previously known p^ . 
Eq j3.39 shows that the cancellation in the denominator in general depends on the relative 
phases of b, c. In addition the solar angle will receive corrections from the charged lepton 
sector from Eq. 2.1C . 

For the case of a partial neutrino mass hierarchy mi ^ 1712 <C arising from the 
right-handed neutrinos having the texture in Eq.|1.5| with the right-handed neutrino with 
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mass Y dominating and the off-diagonal riglit-handed neutrinos of pseudo-Dirac mass X 
giving the sub-dominant contributions, the atmospheric angle is determined similarly to 
the previous case. However in this case the large solar angle arises from a sum of two angles 
as in Eq.3.65, where the angle 9'^2^^'^ is given by Eq. p.66 while the angle is given by 

'12 



Eq,3.67. In this case a large solar angle can result from one or both of the angles 0'^^^^'^ 



9^2 being large. As before cancellations can occur in the denominators, which in general 
depends on the relative phases of 6, c and 6', c' and in addition the solar angle will receive 
corrections from the charged lepton sector from Eq. 2.1C . The pseudo-Dirac limit is not 
relevant for the LMA MSW solution. For both the full and partial hierarchy cases we 
obtain a limit \Ue^\ ~ 0.1, just below the current CHOOZ limit. 

For the case of an inverted neutrino mass hierarchy with a heavier approximately 
pseudo-Dirac neutrino pair, —mi ~ m2 3> rn-3, arising from the right-handed neutrinos 
having the texture in Eq.1.5, with now the off-diagonal pseudo-Dirac right-handed neu- 
trinos of mass X dominating, the large atmospheric angle arises from Eq,3.85, assuming 



the condition in Eq.3.81, similar to ^6[]. In this case the contribution to the solar angle 
from the neutrino sector is almost maximal and we have the result in Eq. 3.104. However, 
including the effect of charged lepton mixing angles, for example a Cabibbo-like contribu- 
tion to 9^2 ~ 0.2, can result in 9i2 in the LMA MSW range and ^13 close to the current 
CHOOZ limit, as in the hierarchical case. 

Although we only quote the leading order expressions, they should serve to provide a 
useful guide to unified model building. More accurate expressions including the 7712/771-3 
corrections can readily be obtained if required from our more general results. 

We have also considered two physical applications of our results, to neutrinoless double 
beta decay and to leptogenesis. Both hierarchical and inverted hierarchical cases predict 
small /?/?oi/ with |777ee| ~ 0.007 eV within the sensitivity of future proposals such as GENIUS. 
We also considered leptogenesis for the sequential sub-dominance case and emphasised that 
successful leptogenesis is possible if the dominant right-handed neutrino is the heaviest one, 
but the leptogenesis phase is unrelated to the MNS phases. Although the leptogenesis phase 
is not directly related to the MNS phases, it should be emphasised that they all originate 
from the same source, namely the Dirac mass matrix. It is natural to assume that the 
neutrino Dirac mass matrix is complex, in which case we would expect both a non-zero 
leptogenesis phase and non-zero MNS phases. The observation of the Dirac MNS phase 
at a neutrino factory would therefore be at least circumstantial evidence for a non-zero 
leptogenesis phase. 

To conclude, recent SNO results when combined with other solar and atmosheric neu- 
trino data supports the LMA MSW solution and three active neutrinos. The application of 
the notion of technical naturalness then leads to two possible mass patterns for neutrinos, 
either a hierarchy or an inverted hierarchy with an approximately pseudo-Dirac heavier 
neutrino pair. Further naturalness requirements applied to the see-saw mechanism, leads 
us to suppose that the successful mass patterns result from right-handed neutrino domi- 
nance where one or two right-handed neutrinos dominate and give the leading order mass 
matrices in Eqs.1.1,1.2. The maximal atmospheric angle results from a ratio of domi- 
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nant right-handed neutrino Dirac couphngs. The large solar angle results from a ratio of 
subdominant right-handed neutrino couplings. The leading order results we present here 
should provide a useful and reliable guide in constructing unified models. The usefulness 
is due to the simplicity of the leading order results and the reliability is due to the correct 
treatment of phases and also the fact that the radiative corrections are expected to be 
relatively small due to the technical naturalness. Both these features are a consequence of 
right-handed neutrino dominance, which may be tested experimentally via its predictions 
of ^13 ^0.1 and \mee\ ~ 0.007, assuming the LMA MSW solution to be correct, although 
we emphasise that both predictions are rather soft since they rely on the assumption of 
the absence of cancellations. 
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A. Equivalence of different parametrisations 



In this appendix we exhibit the equivalence of different parametrisations of the MNS matrix. 
A 3 X 3 unitary matrix may be parametrised by 3 angles and 6 phases. We shall find it 
convenient to parametrise a unitary matrix by ^: 



W = P2R23R13P1R12P3 



(A.l) 



where Rij are a sequence of real rotations corresponding to the Euler angles 6ij, and Pi are 
diagonal phase matrices. The Euler matrices are given by 



where c 



cos 9ij and Sij 



smi 



/l \ 

R23 = C23 S23 
\ -S23 C23 / 

^ Ci3 Si3 

Ri3= 10 

Y -Si3 Ci3 

^ C12 S12 

R12 = -Si2 C12 

V 1 
kj- The phase matrices are given by 

/l o\ 



e'^ 
1 



(A.2) 



(A.3) 



(A.4) 



(A.5) 



^It is convenient to define the paramotrisation of rather than V because the MNS matrix involves 
1/"^^ and the neutrino mixing angles will play a central role. 
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P2 



(l 

e*'^2 












ft 











(A.6) 



(A.7) 



Y 

By commuting the phase matrices to the left, it is not difficult to show that the 



parametrisation in Eq.A.l is equivalent to 

yt = PU23UM2 



where P = P1P2P3 and 





(1 





U23 = 





C23 S23e-'^^-' 






-S23e*'^23 




( 


CIS si3e-*^i3 


Ui3 = 




1 




[- 


5136^^" ci3 




( 


ci2 si2e-*'^i2 


Ui2 = 




si2e'^'^ ci2 




[ 


1 



where 



(^23 = X + ^2 - ^3 
<Jl3 =UJi-UJ3 
612 = UJi- UJ2 



(A.8) 
(A.9) 
(A.IO) 

(A.ll) 

(A.12) 
(A.13) 
(A.14) 



The matrix Umns is an example of a unitary matrix, and as such it may be parametrised 

then the 



by either of the equivalent forms in Eqs. [A.l| or |A.8| . If we use the form in Eq.|i 
phase matrix P on the left may always be removed by an additional charged lepton phase 
rotation AV^'^ = P^ , which is always possible since right-handed charged lepton phase 
rotations can always make the charged lepton masses real. Therefore Umns can always be 
parametrised by 

Umns = U23U13U12 (A.15) 

which involves just three irremoveable physical phases 5ij. In this parametrisation the 
Dirac phase 5 which enters the CP odd part of neutrino oscillation probabilities is given 

by 

(5 = 5^3 - ^23 - 5l2- (A. 16) 



Another common parametrisation of the MNS matrix is 

Umns = R23U13R12P0 



(A.17) 
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where 



( e^^i \ 
1 



(A.18) 



and in Eq.A.17 C/13 is of the form in Eg, A.10| but with (5i3 replaced by the Dirac phase b. 



The parametrisation in Eq. [A.17| can be transformed into the parametrisation in Eq. |A.15 



by commuting the phase matrix Pq iii Eg, A. 17 to the left, and then removing the phases on 
the left-hand side by charged lepton phase rotations. The two parametrisations are then 
related by the phase relations 

(523 = h (A. 19) 

<5i3 = 5 + /?! (A.20) 
5i2 = A - /?2 (A.21) 



The use of the parametrisation in Eg.A.17 is widespread in the literature, however for the 



reasons discussed in the next sub-section we prefer to use the parametrisation in Eg. A. 15 
which is trivially related to Eg.A.17 by the above phase relations. 



B. Charged lepton contributions to the MNS matrix 



In this appendix we discuss the contribution of the charged lepton mixing angles to the 
MNS matrix. The MNS matrix is constructed in Eg.^j^ as a product of a unitary matrix 
from the charged lepton sector and a unitary matrix from the neutrino sector V^^^ . 
Each of these unitary matrices may be parametrised by the parametrisation of V'^ in Eg.A.l 
Thus we write 



- p'^L rft-'L n'^L r>>^L n'^L id^l 

- ^2 -"^23 -"-13 ^1 ^12 ^3 

pEL pEi^ P^L pE^L p-^L p-El 
^2 ^23 -"-13 -^1 -"-12 -^3 



(B.l) 
(B.2) 



where the Euler angles and phases are defined as in Egs. A.2-A.7 but now there are indepen- 
dent angles and phases for the left-handed neutrino and charged lepton sectors distinguished 
by the superscripts and E^. The MNS matrix from Egs.p^ B.l, p.2 is then 



TJ _ pEL^ pEL^ pEL^ pEL^ pEL^ pEL^ p^L pl^L pl^L pl^L p^L p 

UmNS - ^3 ^12 ^1 ^13 ^23 ^2 ^2 ^23 ^13 ^1 ^12 ^ 



3 



(B.3) 



E t 

As before we commute all the phase matrices to the left, then choose Pg ^ to cancel all 
the phases on the left-hand side, to leave just 



(B.4) 



with independent phases and angles for the left-handed neutrino and charged lepton sectors. 



in the convention of Egs . |A . 9| , |A . 1 0| , |A . 1 1| . The phases in Eg. |B.4| are given in terms of the 



phases in Egs. B.l, B.2 by 



5S 



UJ. 



UJ 



VL 



UJ 



VL 



(B.5) 
(B.6) 
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"13 
"12 



^2^ 



■ IjJ. 



(B.7) 
(B.8) 
(B.9) 

X- + 02 " - '/'2" - X-"' + - (B-10) 
The form of Umns in Eq jB.4 is similar to the parametrisation in Eq. 2^, which is the 



practical reason why we prefer that form rather than that in Eq. A.17 . 

We now discuss the MNS matrix to leading order in ^13. From Eqs. A.15, A.£. AlTo|jA.ll| , 
we find to leading order in ^13 that Umns may be expanded as: 



Umns 
/ 

\ Si2S23e 



C12 



Sue 



Sl2C23e 

i{S23+Si2) 



Sl2S23013e^(-^23+5i3-5i2) 5236-^^23 



Cl2C23D'l3e 



-Cl2S23e 



iS23 



Sl2C23P13e 



j((5l3— <5l2) 



C23 



(B.ll) 



For ^13 = 0.1, close to the CHOOZ limit, the approximate form in Eq. |B.ll| is accurate to 
1%. 

We now wish to expand the MNS matrix in terms of neutrino and charged lepton 



mixing angles and phases to leading order in small angles, using Eq.B^. In technically 
natural theories, based on right-handed neutrino dominance, the contribution to ^23 comes 
mainly from the neutrino sector, ^23 ~ ^23 • Furthermore in natural theories we expect that 
the contributions to ^13 are all separately small so that the smallness of this angle does 
not rely on accidental cancellations. Clearly this implies that 9'^^ and 9^^ must both be 
^ 013 . Since the 13 element of Umns also receives a contribution from the charged lepton 
sector proportional to s^^Sgs; the same argument also implies that 9^2 ~ ^13- Therefore 
the natural expectation is that all the charged lepton mixing angles are small! Expanding 



Eq jB.4| to leading order in small angles 9^2 1 ^23 ' ^13 ' ^13 ' 



Umns 



-12 



'12' 



™12 



'"13 



-^23) r^'.kr'ii^ - .^.k.^^kft':ip^^-^23+Q~^'i2) .^kp-^^23 



•^12*23^ 



^i^L „<^L a<^L piS'^k 
•-12 '-23 ^^13 ^ 



'12 "-23 

-C^kSoke^^23 



:)El 
'23 



+9 



El 
13 



+^12 



-12 



'12^ 



'12 •'23' 

™12 



'23 ''12' 
^23^12*^ 



-23 



^°23 



+5 2^) 



-12 •'23' 



'12 



-23 •'12 

V 

"12 



'12 



-23 ^-12 



J^L „-ih 
'^23^ 










'23" 



«12C23^13e^(''l3 "^12 ) 



-23 
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(B.12) 



where we have dropped terms of order 62^013. The first matrix on the right hand side of 



Eq jB.12| gives the contribution to the MNS matrix from the neutrino mixing angles and 
phases, and is of the same form as Eq jB.ll . The subsequent matrices give the corrections 
to the MNS matrix from the charged lepton mixing angles 623 , O13 , and 6^2 ■ 



C. Proceedure for diagonalising hierarchical mass matrices 

In this appendix we discuss the diagonalisation of a general complex hierarchical matrix 
m where 

^ mil mi2 mi3 
m = I m2i m22 "^23 (CI) 
msi m32 mss 

The matrix m is diagonalised by a sequence of tranformations: 



R 



L* r>L 



^12 



T 



T 



nij " r>iy " r>iy -pR jyR jyR pR pR pR 
^13 ^23 ^2 "1-^2 ^23^13-^1 ^12^c 



I mi Q 
m2 
ma 



(C.2) 



In the case of the charged lepton mass matrix, all the rotation angles are small, while in 
the case of the neutrino mass matrix it is symmetric. The results for the general complex 
matrix m will be sufficiently general to allow us to apply them to both of the physical cases 
of interest as limiting cases. 

The proceedure for diagonalising a general hierarchical matrix m involves the following 
steps. 

1. The first step involves multiplying the mass matrix m by the inner phase matrices 



P2 defined in Eg. A. 6 
/ 



P^'mPi 



R 



mil 



m2ie 
\m31e- 



mi2e^2 



m22e- 
m32e* 



miae 

i 



m23e 



mil "i'i2 "i'i3 



m2i 
m'q 



m22 m' 



23 



(C.3) 



The purpose of this re-phasing is to facilitate steps 2,3 using real rotation angles 623, ^13, 
as we shall see. 



2. The second step is to perform the real rotations ii23 defined in Eq jA.2| on the 
re-phased matrix from stepl. The purpose is to put zeroes in the 23,32 elements of the 
resulting matrix: 



^23 



( 



mil f^i2 "i'i3 



m2i m. 



22 ^23 



m^i m 



32 ^33 



R23 



( mil mu rhis 

77121 "122 



(C.4) 



\"131 







m. 



The zeroes in the 23,32 positions are achieved by diagonalising the lower 23 block, using 
the reduced matrix i223 obtained by striking out the row and column in which the unit 



- 31 - 



element appears, to leave a 2 x 2 rotation, 



^23 



17122 



23 



^23 



m22 
m'o 



(C.5) 



which implies 



tan 20^ 



2 [m^jgm'aa + m^2'^32] 



23 



tan26'|| 



m' 



I 2 

33 
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I 2 

22 



+ m 



I 2 

32 



m 



I 2 

23 



2 [m'33m^j2 + m'aam'ag] 



m' 



I 2 

33 



/ 2 

22 



+ m: 



/ 2 

23 



/ 2 

32 



(C.6) 



(C.7) 



The requirement that the angles ^23 ^23 ^^^^ means that the numerators and 
denominators must have equal phases, and this is achieved by adjusting the relative phases 
(f)f- — (pj which appear in the lower block of Eq. p.3| . The remaining elements are then given 
by the reduced rotations 

(C.8) 
(C.9) 



mi2 mi3 




3. The third step is to perform the real small angle rotations R13 defined in Eg, A.3 on 
the matrix from step 2. The purpose is to put zeroes in the 13,31 elements of the resulting 
matrix: 

mil rhi2 

77121 ?7l22 



R 



13 



\rri31 



m 



pi? 

^13 



3/ 



m21 77122 

\ m'^ 



(C.IO) 



The zeroes in the 13,31 positions are achieved by diagonalising the outer 13 block, using 
the reduced matrix Ri^ obtained by striking out the row and column in which the unit 
element appears, to leave a 2 x 2 rotation. 



R 



L T I mil mi3 
\ 77731 7773 



pi? 

^13 



mil 
m'n 



(C.ll) 



which implies 



'13 



77713 ^-317^11 



7773 

777n 



+ 



(777^5)2 
,f\2 



(C.12) 



(C.13) 



.3 (777^)2 

The requirement that the angles 6^^ and 9^^ are real fixes the absolute value of the phases 
+ since only the relative phases were fixed previously. This uses up all the phase 
freedom and thus all the resulting mass matrix elements in Eq JC.lC remain complex. Note 
that Eq. p.lC is written to leading order in the small angles 9is , and as discussed previously 
the 23,32 elements remain zero to this order. The large complex element 7^3 is approx- 
imately unchanged to this order. Due to the zeroes in the 23,32 position of the matrix 
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the elements mi2 and 77121 are also unchanged to leading order. The element 77122 is also 
unchanged of course since it is not present in the reduced matrix. The only new element 
is therefore 

7711377131 



7^11 W 77711 



777q 



(C.14) 



4. The fourth step involves multiplying the mass matrix resulting from Eq. C.10| by the 



phase matrices Pi defined in Eq. |A.5 
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The purpose of this re-phasing is to facilitate step 5 using real rotation angle ^12. 



5. The fifth step is to perform the real rotations P12 defined in Eq. on the re-phased 
matrix from step 4. The purpose is to put zeroes in the 12,21 elements of the resulting 
matrix: 
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The zeroes in the 12,21 positions are achieved by diagonalising the upper 12 block, using 
the reduced matrix P12 obtained by striking out the row and column in which the unit 
element appears, to leave a 2 x 2 rotation, 
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(C.17) 



which implies 



tan 26'f. 
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(C.19) 



'-21 y 



The requirement that the angles 6^2 and 621 are real means that the numerators and 
denominators must have equal phases, and this is achieved by adjusting the phases xl, Xr- 
6. The sixth step involves multiplying the complex diagonal mass matrix resulting 
from Eq. C.16| by the phase matrices P3 defined in Eq.| 
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The result of this re-phasing is a diagonal matrix with real eigenvalues. In the case of 
charged leptons this last step can be achieved by a suitable for any choice of P^. This 
freedom in P^ enables three phases to be removed from the MNS matrix. 
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D. Diagonalising the hierarchical neutrino mass matrix 



In this appendix we shall apply the results of appendix ^ to the case of the complex 
symmetric hierarchical neutrino mass matrix of the leading order form as shown in Eq. |l.l| , 
which will be written in full generality as 
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LL 
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\m 
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12 1 
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13 
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12 |mVq|e*"'^i3 



12 m: 



221 



e*<^i3 |m23le'"^23 \m 




(D.l) 



where it should be remembered that the elements in the lower 23 block are larger than the 



other elements, as in Eq.1.1 



The proceedure outlined in appendix ^ for diagonalising ?ti^^ is to work our way from 
the inner transformations to the outer transformations as follows. 



1. Re-phase m^^;^ using the Pr 

2. Put zeroes in the 23=32 positions using -R23- 

3. Put zeroes in the 13=31 positions using R^k,- 



4. Re-phase the mass matrix using 



5. Put zeroes in the 12=21 positions using 

6. Make the diagonal elements real using the Pi 



3 • 

If is small, then for the hierarchical case mj, ^ m2 this proceedure will result in an 
approximately diagonal matrix to leading order in ^^3 . One might object that after step 
3 the i2^3 rotations will "fill-in" the zeroes in the 23,32 positions with terms of order Q^^ 



multiplied by m 



12 ' '"'13 • 



However in the hierarchical case m 



12 ' "^13 



are smaller than m' 



13 

33 



by a factor of 9-^^ which means that the "filled-in" 23,32 entries are suppressed by a total 
factor of (^13)^ compared to the 33 element. This means that after the 5 steps above a 
hierarchical matrix will be diagonal to leading order in ^^3 , as claimed. For the inverted 
hierarchical neutrino case a different proceedure must be followed, as discussed in the next 
sub-section. Here we shall systematically diagonalise the hierarchical neutrino mass matrix 



in Eq.D.l by following the above proceedure as follows. 



The first step is to re-phase the matrix in Eq. D.l using P^^* so that the neutrino mass 
matrix becomes, 
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(D.2) 



To determine the 23 neutrino mixing angle ^23 we perform a 23 rotation which diagonalises 



the lower 23 block of Eq. p.2| . From Eq.C.C we find the 23 neutrino mixing angle ^23 



The relative phase 
real, 



is fixed by the requirement that the angle ^23 Eq.D.3 be 
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After the 23 rotation in Eq.C.4, the neutrino mass matrix in Eq. D.2 becomes 
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The lower block elements are given by 
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which implies 
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We now perform a 13 rotation on the neutrino matrix in Eq.D.5 which diagonalises the 



outer 13 block of Eq.D.5 and determines the 13 neutrino mixing angle d^k- From Eq.C.12 



we find the 13 neutrino mixing angle as 
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The absolute phases 
be real, 



62^, are fixed by the requirement that the angle in Eq. D.lO 
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After the 13 rotation in Eq.C.lO, Eq.p.5| becomes 
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To leading order in the only new element in Eq. p.l2| is 
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It only remains to determine the 12 neutrino mixing angle by diagonalising the 



upper 12 block of Eq,D.12. From Eq.C.lS we find the 12 neutrino mixing angle as 
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The phase x"^^ is fixed by the requirement that the angle 9'^2 Eq.D.14 be real, 
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After the 12 rotation the upper block of the matrix in Eq. D.12 is diagonal and the resulting 
matrix is 

m[ \ ( mie^'t>'i 

m2e^'t''2 I (D.16) 
m3e^<^3 



m'2 
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where from Eq. p.l7 



m'2 = (s5'f)Vn|e''^" +2sSc5't|mi2|e'('^"-^''^) + (c^'^) V22|e'^^^^"'''"'^^ 



(D.17) 
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It is a simple matter to adjust the phases uj^^ in Pg^ to remove the phases in Eq. p.l6| 

(D.19) 



and make the neutrino masses real, as in Eq.C.20 



I o 



This completes the diagonalisation in Eq.3.3. In the case of neutrino masses, unlike the 
case of the charged fermions, there is no left over phase freedom. This is the reason why 
the MNS matrix has three more physical phases than the CKM matrix. 

E. Diagonalising the inverted hierarchical neutrino mass matrix 

In this appendix we shall consider the case of the complex symmetric inverted hierarchical 
neutrino mass matrix of the leading order form as shown in Eq.^]^. In this case the 
proceedure is as follows. 

1. Re-phase m'^]^ using the P2^. 

2. Put zeroes in the 13=31 positions using i?23- 

3. Put zeroes in the 23=32 positions using Ri^- 

4. Re-phase the mass matrix using P^^ . 

5. Put zeroes in the 12=21 positions using Ri2- 
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6. Make the diagonal elements real using the Pi 



We continue to write the neutrino mass matrix as in Eq.D.l, but now it should be 



remembered that the 12,13 elements are now larger than the other elements, as in Eq.|L2 



This is the reason why the above proceedure differs from that for the case of the hierarchical 
neutrino mass matrix. 



We first perform the re-phasing as in Eq.D.2. Then we determine the 23 neutrino 
mixing angle by performing a 23 rotation such that 
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which gives the 23 neutrino mixing angle ^23 this case to be 
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Since the Euler angles are constrained to satisfy 
and this then fixes 



< 7r/2, we must have tan ^23 ~ +1 
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and gives 
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The lower block elements are given by 
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which implies 



(|m^2|e'^'^^^"'^^^^ - |m^3|e'(<^33-203^)) + ((c-^)2 _ (^f )2)|m^3|e'(^23-<^2' 



and the remaining diagonal elements are given as before in Eqs.D.7,D. 
We next perform a small angle 13 rotation such that 
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where 
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Note that to leading order in 6\i^ the large element is unchanged. The remaining 
elements in Eq.^ 

which is of order 



are also unchanged to leading order in O'^i^. The 13=31 element in 



Eq,E.£ gets filled in by a term 6 



m 



11 



^13)^ compared to rh'^ 



12 



and does not appear to leading order in From Eq. E.lO] the 13 neutrino mixing angle 
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The requirement that 9^!^ is real fixes the absolute value of the phases 02 



The left hand side of Eq.E.E now resembles the left hand side of Eq.D.12, except that 
here m^^ is unchanged due to the zero 13=31 element after the 23 rotation. Therefore 
the rest of the diagonalisation process follows that of the previous hierarchical case from 



Eq,D.14 onwards, where now 



tan29if 
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Note that in the inverted hierarchy case here we have 

I"il2l > l"^22l> l"^lll (E.13) 

which implies an almost degenerate pair of pseudo-Dirac neutrino masses (with opposite 



sign eigenvalues) and an almost maximal 12 mixing angle from Eq. E.12 
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